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GROUND  STATE  OF  ATOMIC  HELIUM 

ABSTRACT 

A  variational  calculation  of  the  ground  state  energy  Is 
made  using  basis  functions  with  the  correct  behavior  at  the  simple 
two-electron  singularities.   The  truncation  errors  obtained  with 
few  terms  are  not  orders  of  magnitudes  less  than  those  obtained 
with  a  spherical  harmonics  expansion.   It  seems  that  the  slow 
convergence  of  the  Legendre  expansion  is  not  due  to  its  inability 
to  account  for  this  singularity  with  only  a  finite  number  of  terms 
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GROUND  STATE  OF  ATOMIC  HELIUM 

1 .    Introduction. 

There  is  a  long  history  of  Rayleigh-Ritz  calculations  of  the 
lowest  energy  eigenvalue  of  the  helium  atom  and  isoelectronic 
ions.    Many  contemporary  studies  have  special  purposes;   one  is 
to  test  trial  functions  with  counterparts  in  other  problems. 
In  particular,  we  want  to  avoid  trial  functions  with  a  complicated 
dependence  upon  r  ^      (the  interelectronic  distance),  which  are 
used  in  high-accuracy  calculations  [2,5,^]-   One  attempt  is  to 
avoid  the  use  of   r,p   completely;  this  often  leads  to  a  Legendre 
expansion  [5].   At  first  it  seemed  to  converge  rapidly,  but  later 
studies  indicate  a  more  complicated  behavior  [6,7,8,9].   The  first 
few  terms  do  indeed  converge  rapidly,  but  they  converge  to  a  value 
bounded  from  the  true  eigenvalue.   After  this  the  expansion  converges 
so  slowly  it  is  useless. 

This  unexpected  behavior  of  the  Legendre  expansion  is  not 
completely  understood.   It  has  been  conjectured  that  this  is  related 
to  the  singularities  of  the  Schroedlnger  equation.   Specifically, 
the  continuity  properties  of  the  derivatives  of  the  actual  Schroedlnger 
function  are  different  from  those  of  the  individual  members  of  its 


For  an  excellent  summary,  see  Bethe  and  Salpeter  [l] 

-  h  - 


Legendre  expansion.   Consequently,  the  derived  series  must  not 

converge  uniformly,  and  the  coefficient  of  the  P. —  Legendre 

-2  * 
polynomial  cannot  decrease  faster  than   ^   .    We  wanted  to  test 

the  conjecture  by  performing  a  Raylelgh-Rltz  calculation  using  a 

trial  function  constructed  so  that  a  finite  (small)  number  of 

terms  might  have  the  same  behavior  at  the  singularities  as  the 

solution  of  the  Schroedlnger  equation.   We  could  not  do  this 

completely;  the  singularity  due  to  the  coincidence  of  all  three 

charges  Is  too  complicated.   However  It  seemed  worthwhile  to 

do  a  more  restricted  calculation  which  takes  account  of  the 

simple  singularities. 

At  present  we  have  only  preliminary  results,  which  are 

essentially  negative.   Our  best  energy  eigenvalue   approximately 

equals  the   F-llmlt,  I.e.  the  best  eigenvalue  obtainable  with  a 

Legendre  expansion  Including  spherical  harmonics  of  orders  no 

greater  than  three  [7  ] .   This  result  Indicates  that  either  the 

large  scale  Interelectron  correlations  cannot  be  represented  by 

simple  wave  functions,  or  the  complicated  singularities  due  to 

coincidence  of  more  than  two  charged  particles  Is  about  as 

Important  as  the  simple  ones.   The  first  alternative  Is  somewhat 

unwelcome,  for  although  the  error  In  the  energy  can  easily  be 

reduced  to  less  than   10"-^   times  the  Ionization  potential  for 

atomic  helium,  large  scale  correlations  are  expected  to  be  more 

Important  In  many  other  problems . 


Schwartz  '[8]  obtained  the  corresponding  result  for  the  energy, 

-  R  - 


II.   General  Theory. 

The  Schroedlnger  equation  for  a  two-electron  atom  Is 
2 


(2.1) 


1 


^-|vM,). 


1=1 


12 


g!  =  E0 


We  use  the  usual  notation;   Z   is  the  nuclear  charge  at  the  origin, 

r.   the  radius  of  the   1 —  electron,  and   r..^   the  interelectronic 
1  12 

distance.   Equation  (2.l)  and  the  condition  0     is  in  Lp   do  not 
define  an  eigenvalue  problem,   Kato  has  shown  there  exists  a  self- 
adjoint  problem  such  that  equation  (2.1)  holds  at  its  ordinary 
points  (where  no  potential  diverges)  [lO],  and   0  is  everywhere 
continuous  [ll].   The  local  properties  of  0     follow  from  (2.1) 
and  continuity  alone.*   When  v-.^      vanishes  and   r.   does  not 
vanish,  the  singularity  is  easily  characterized.   We  Introduce 
as  new  arguments  of  0     vectors   r   and  r-.^,   respectively  the 
electrons'  center-of-mass  and  relative  coordinate. 


(2.2) 


12     2 

T  ^o  -  ^12 


z(^  +  2_)  +  J^ 

^1    ^2     ^12 


0  =   E0 


We  expand  0     formally  in  powers  of   r-,p   and  substitute  into 

(2.2)  and  obtain 

(2.3)  0  =  a^  +   a^r^g  +  ^(r^g^ 


(2.4) 


-  2a,  +  a   =  0 
1    o 


We  are  Indebted  to  Prof.  P.  John  for  discussions  concerning 
this  point. 
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As  a  result  of  (2.4),   0  may  be  written  In  the  form 

(2.5)  0  =  0j  +  I  r^2  ^11 

where  0j     and  0^^  are  equal  if  r,p  vanishes,  and  are  analytic* 
in  any  domain  In  which  neither  r.   vanishes. 

If  00  denotes  the  cosine  of  the  angle  "between  the  radius 
vectors  of  the  electrons,  the  usual  Legendre  expansion  L5j  is 

oo 

(2.6)  0  =  \~     P^Cr^r2)P^(aO 

£=0 
The  functions  0j     and  0^.^  may  be  expanded  similarly.   The 
relation  between  these  series  involves 

£  J+2     \ 

CO    /    T      r  .         1  ^<  \ 

^^   ^  ^s^  2/.-1  vj    ^        2£+J>   r^^^"-  /   ^ 

Here  r^  and  r^   denote  respectively  the  lesser  and  greater  of 
r.,   and  r^ ,   If  the  series  for  0j     and  0^^  converge  rapidly, 
the  expansion  for   0  must  converge  slowly.   The  crucial  point 
is  that  when   r.p   vanishes,  the  analytic  properties  of  the 
individual  terms  of  (2,6)  differ   from  those  of   0  itself.   In 
contrast  to  this,  the  individual  terms  of  (2,6)  have  the  correct 
behavior  when  one  electron  is  at  the  origln>  this  singularity  does 
not  affect  the  convergence , 

We  construct  a  trial  function  0     to  approximate  0     following 
(2.5) 

A  function  is  called  anytlc  if  it  is  an  analytic  function  of  the 
six  Cartesian  coordinates  of  the  individual  electrons. 
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(2.8)      ?  =  ?j  +  I  r^2^II 


(2.9a)     ?T  =  V~  RT^fri^QHr,r„)^P,(a)) 


Ii'"l  2"  1  2'  H 

I 

(2.9b)     \^   =   ^  Rm^ri^2^^^1^2-^^Pi'^^)    • 

It  Is  convenient  to  represent  the  radial  functions  by  exponentials, 

(2.10a)    R;^^  =  (l+P^g)  XI^I^v^^P^->^iv^l-^2^v^2^ 

V 

(2.10b)     R^^j^   =  {1+?^^)   ^ajj^^^expi-K^^^r^-K^^.g^g) 

V 

Here   P-,p   Is  an  operator  which  Interchanges   r,   and  r^,   and 
a's   and  M's   with  various  subscripts  are  variational  parameters. 
We  restrict   ^  to  satisfy  (2,4)  as  follows:  , 


(2.11)     a^j^g^  -  a-j.^^ 


(2.12)     K^^,  +K2^^  ^X^^,  +V<^^^ 


These  restrictions  are  more  stringent  than  needed,  but  conveniently 
give  an  expansion  whose  Individual  terms  satisfy  (2.4). 
The  trial  energy  Is 


(2.15)  E   = 


■^  ir?TJi9?^°P 


%2,6 


y^d 


r 


Here   H   Is  the  operator  on  the  left-hand  side  of  (2.1)  and  the 
Integration  variables  are  the  Cartesian  coordinates  of  the 
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Individual  electrons.   The  formulae  for  all  integrals  needed  to 
evaluate  (2.13)  are  routine;  the  only  one  that  might  at  first 
seem  otherwise  is 

An  advantage  of  radial  functions  of  the  form  (2.10)  is  that  all 
integrals  are  the  same  functions  of  different  arguments.   We  do 
not  write  out  an  explicit  formula  for  the  right-hand  sides  of 
(2.13)  as  it  is  too  complicated  to  be  illuminating. 

III.  Numerical  Results  and  Discussion. 

Using  the  IBM-709O  electronic  digital  computer,  we  evaluated' 
E  for  atomic  helium  with  different  sets  of  values  of  the  variational 
parameters.   The  values  of  all  integrals  involved  were  carefully 
checked  by  comparing  the  machine  results  for  special  cases  with 
claculations  using  a  desk  computer.   The  formulae  for  these  cases 
were  derived  independently.   We  believe  it  is  very  unlikely  that 
our  values  of   E   for  a  particular  set  of  variational  parameters 
are  incorrect.   However  the  problem  of  programming  the  computer 
to  search  for  the  minimum  of  E  proved  complicated.   We  reached 
the  stage  where  a  great  deal  of  additional  work  is  needed  to 
improve  this  part  of  the  program  significantly.*   On  the  other 
hand,  the  partial  results  already  obtained  warrant  communication. 


The  difficulty  of  minimizing  programs  is  that  modern  computers 
are  not  suited  to  making  complex  decisions.   If  a  well  understood 
problem  like   He   is  difficult,  we  must  doubt  results  in  problems 
not  well  understood. 
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The  accuracy  already  obtained  may  be  Interesting  since  the  trial 
functions  have  relatively  few  terms.   In  view  of  the  stationary 
nature  of  E,   we  believe  that  the  best  values  of  the  parameters 
will  not  reduce  the  errors  by  orders  of  magnitude. 

If  only  terms  with  i   =   0     are  retained  In  (2.9),  we  have 
results  for  three  terms  in  '2.10).   We  get 

(5.1)  E  =  -  2.89781 

The  difference  between  this  and  our  best  value  for  only  two  terms 
in  (2.10)  is  only   O.OOO3O,   which  suggests  that  increasing  the 
number  of  terms  in  (2.10)  will  not  lead  to  much  lower  values  of 
E.   The  above  value  may  be  compared  to  the  S-llmit  [7j8] 

(3.2)  E(0)  =  -  2.87902 

This  comparison  is  not  illuminating  because  the   r,^   terms  in 
our  trial  function  contain  considerable  angular  anti-correlation; 
this  rather  than  the  effect  of  the  singularity  probably  accounts 
for  the  difference.   If  terms  with  £   -^  1      are  retained  in  (2.9), 
we  get  :•...■;;-.. 

(5.3)  E  =  -  2.90313     ■ •   . 

This  may  be  compared  to  estimated  limits  if  terms  with   -^  ;1  1 
or   ^  <_  3   are  retained  in  [2.6)    [7] 

(3.^)  E(l)  =  -  (2.90049  +  0.0001) 

(3.5)  E(3)  =  -  (2.90332  +  0.0001) 

Despite  some  uncertainty  in  the  estimated  limits   the  value  of 
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E  as  given  by  (3 -5)  is  disappointing.   If  the  value  of  E 
cannot  be  Improved  without  going  to  very  high   i   In  (2.10)  we 
are  led  to  the  two  alternatives  mentioned  In  the  Introduction. 
Either  the  large-scale  electron  correlations  give  rise  to  the 
dominant  errors  In  (:5'3)  and  (3-5)>  or  the  complicated  singularities, 
which  arise  from  coincidence  of  more  than  two  charges,  are  as 
important  In  their  effects  upon  E  as  the  simple  ones.   The  second 
alternative  Is  hard  to  believe.   If  the  first  alternative  Is  correct, 
we  expect  that  the  differences  between  the  exact*  value  of  E   and 
the  limiting  values   E(^)   would  continue  to  decrease  rapidly  when 
i  ^  5)   in  contradiction  to  the  results  of  Schwartz  [8]. 

We  want  to  note  a  few  features  of  our  trial  function   ^. 
Suppose  we  keep  only  the  one  term  In  (2,9)^  (2.10). 

(3.6)  V  -    (I+Pi2)(e   112  2  +  1  p^^e   5  1   ^  2^ 

(3.7)  Ki  +  V<2  =  H3  +  X4 


The  optimum  parameters  and  energies  are 
(3.8)  E  =  -  2.896665 

(3.9a)  M^  =   2..351 

r 

(3.9b)  K2  =  i-374 

(3.9c}  H^  =  1.927 

(3.9d)  M^  =  1.798 


*  For  our  purposes,  we  assume  this  Is   -  2.903724,  the  six  decimal 
approximation  to  Pekerls '  result  [4], 
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The  eigenvalue  (3>8)  Is  good  for  a  simple  wa^z-e  function.   Now 
consider  the  asymptotic  case   r^  ^^  r. 


^         -X^v^    -  H^v, 


(5.10)  0  ^e  ^  ^    ^  2 
The  correct  asymptotic  form  is 

(5.11)  0-^r2~°-^^^  exp(-2r^-1.5^^r2) 

The  most  striking  failure  of  (5,10)  Is  that  the  Inner  electron 
Is  too  close  to  the  origin.   Finally,  we  note  that  our  best  trial 
function  thus  far,  with  five  terms  In  (2.10)  whose  parameters 
are  listed  In  table  I,  gives  a  much  worse  asymptotic  formula 

^     -l,808r,  -  0.525ro 

(5.12)  0-e        1         2 

Perhaps  this  indicates  how  little  the  asymptotic  form  affects   E. 
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TABLE   I 

s, 

V 

Xl 

Kg 

>^3 

^4 

a 

0 

1 

2.351 

1.374 

1.927 

1.798 

7.474 

0 

2 

3.609 

3.609 

3.609 

3.609 

1.355 

0 

5 

1.808 

0.523 

1.165 

1.165 

0.365 

1 

1 

2.515 

2.515 

2.515 

2.515 

2.255 

1 

2 

5.91^ 

5.914 

5»9l4 

5»9i4 

7.799 
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